3.3.36 \(\int \frac {x^2 (a+b \text {arcsinh}(c x))^2}{(d+c^2 d x^2)^2} \, dx\) [236]

3.3.36.1 Optimal result
3.3.36.2 Mathematica [A] (verified)
3.3.36.3 Rubi [A] (verified)
3.3.36.4 Maple [F]
3.3.36.5 Fricas [F]
3.3.36.6 Sympy [F]
3.3.36.7 Maxima [F]
3.3.36.8 Giac [F]
3.3.36.9 Mupad [F(-1)]

3.3.36.1 Optimal result

Integrand size = 26, antiderivative size = 213 \[ \int \frac {x^2 (a+b \text {arcsinh}(c x))^2}{\left (d+c^2 d x^2\right )^2} \, dx=-\frac {b (a+b \text {arcsinh}(c x))}{c^3 d^2 \sqrt {1+c^2 x^2}}-\frac {x (a+b \text {arcsinh}(c x))^2}{2 c^2 d^2 \left (1+c^2 x^2\right )}+\frac {(a+b \text {arcsinh}(c x))^2 \arctan \left (e^{\text {arcsinh}(c x)}\right )}{c^3 d^2}+\frac {b^2 \arctan (c x)}{c^3 d^2}-\frac {i b (a+b \text {arcsinh}(c x)) \operatorname {PolyLog}\left (2,-i e^{\text {arcsinh}(c x)}\right )}{c^3 d^2}+\frac {i b (a+b \text {arcsinh}(c x)) \operatorname {PolyLog}\left (2,i e^{\text {arcsinh}(c x)}\right )}{c^3 d^2}+\frac {i b^2 \operatorname {PolyLog}\left (3,-i e^{\text {arcsinh}(c x)}\right )}{c^3 d^2}-\frac {i b^2 \operatorname {PolyLog}\left (3,i e^{\text {arcsinh}(c x)}\right )}{c^3 d^2} \]

output
-1/2*x*(a+b*arcsinh(c*x))^2/c^2/d^2/(c^2*x^2+1)+(a+b*arcsinh(c*x))^2*arcta 
n(c*x+(c^2*x^2+1)^(1/2))/c^3/d^2+b^2*arctan(c*x)/c^3/d^2-I*b*(a+b*arcsinh( 
c*x))*polylog(2,-I*(c*x+(c^2*x^2+1)^(1/2)))/c^3/d^2+I*b*(a+b*arcsinh(c*x)) 
*polylog(2,I*(c*x+(c^2*x^2+1)^(1/2)))/c^3/d^2+I*b^2*polylog(3,-I*(c*x+(c^2 
*x^2+1)^(1/2)))/c^3/d^2-I*b^2*polylog(3,I*(c*x+(c^2*x^2+1)^(1/2)))/c^3/d^2 
-b*(a+b*arcsinh(c*x))/c^3/d^2/(c^2*x^2+1)^(1/2)
 
3.3.36.2 Mathematica [A] (verified)

Time = 1.47 (sec) , antiderivative size = 385, normalized size of antiderivative = 1.81 \[ \int \frac {x^2 (a+b \text {arcsinh}(c x))^2}{\left (d+c^2 d x^2\right )^2} \, dx=-\frac {\frac {a^2 c x}{1+c^2 x^2}+\frac {2 b^2 \text {arcsinh}(c x)}{\sqrt {1+c^2 x^2}}+\frac {b^2 c x \text {arcsinh}(c x)^2}{1+c^2 x^2}+\frac {a b \left (-i \sqrt {1+c^2 x^2}+\text {arcsinh}(c x)\right )}{-i+c x}+\frac {a b \left (i \sqrt {1+c^2 x^2}+\text {arcsinh}(c x)\right )}{i+c x}-a^2 \arctan (c x)-\frac {1}{2} i a b \left (\text {arcsinh}(c x) \left (\text {arcsinh}(c x)-4 \log \left (1+i e^{\text {arcsinh}(c x)}\right )\right )-4 \operatorname {PolyLog}\left (2,-i e^{\text {arcsinh}(c x)}\right )\right )+\frac {1}{2} i a b \left (\text {arcsinh}(c x) \left (\text {arcsinh}(c x)-4 \log \left (1-i e^{\text {arcsinh}(c x)}\right )\right )-4 \operatorname {PolyLog}\left (2,i e^{\text {arcsinh}(c x)}\right )\right )+i b^2 \left (4 i \arctan \left (\tanh \left (\frac {1}{2} \text {arcsinh}(c x)\right )\right )+\text {arcsinh}(c x)^2 \log \left (1-i e^{-\text {arcsinh}(c x)}\right )-\text {arcsinh}(c x)^2 \log \left (1+i e^{-\text {arcsinh}(c x)}\right )+2 \text {arcsinh}(c x) \operatorname {PolyLog}\left (2,-i e^{-\text {arcsinh}(c x)}\right )-2 \text {arcsinh}(c x) \operatorname {PolyLog}\left (2,i e^{-\text {arcsinh}(c x)}\right )+2 \operatorname {PolyLog}\left (3,-i e^{-\text {arcsinh}(c x)}\right )-2 \operatorname {PolyLog}\left (3,i e^{-\text {arcsinh}(c x)}\right )\right )}{2 c^3 d^2} \]

input
Integrate[(x^2*(a + b*ArcSinh[c*x])^2)/(d + c^2*d*x^2)^2,x]
 
output
-1/2*((a^2*c*x)/(1 + c^2*x^2) + (2*b^2*ArcSinh[c*x])/Sqrt[1 + c^2*x^2] + ( 
b^2*c*x*ArcSinh[c*x]^2)/(1 + c^2*x^2) + (a*b*((-I)*Sqrt[1 + c^2*x^2] + Arc 
Sinh[c*x]))/(-I + c*x) + (a*b*(I*Sqrt[1 + c^2*x^2] + ArcSinh[c*x]))/(I + c 
*x) - a^2*ArcTan[c*x] - (I/2)*a*b*(ArcSinh[c*x]*(ArcSinh[c*x] - 4*Log[1 + 
I*E^ArcSinh[c*x]]) - 4*PolyLog[2, (-I)*E^ArcSinh[c*x]]) + (I/2)*a*b*(ArcSi 
nh[c*x]*(ArcSinh[c*x] - 4*Log[1 - I*E^ArcSinh[c*x]]) - 4*PolyLog[2, I*E^Ar 
cSinh[c*x]]) + I*b^2*((4*I)*ArcTan[Tanh[ArcSinh[c*x]/2]] + ArcSinh[c*x]^2* 
Log[1 - I/E^ArcSinh[c*x]] - ArcSinh[c*x]^2*Log[1 + I/E^ArcSinh[c*x]] + 2*A 
rcSinh[c*x]*PolyLog[2, (-I)/E^ArcSinh[c*x]] - 2*ArcSinh[c*x]*PolyLog[2, I/ 
E^ArcSinh[c*x]] + 2*PolyLog[3, (-I)/E^ArcSinh[c*x]] - 2*PolyLog[3, I/E^Arc 
Sinh[c*x]]))/(c^3*d^2)
 
3.3.36.3 Rubi [A] (verified)

Time = 1.19 (sec) , antiderivative size = 191, normalized size of antiderivative = 0.90, number of steps used = 11, number of rules used = 10, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.385, Rules used = {6225, 27, 6204, 3042, 4668, 3011, 2720, 6213, 216, 7143}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {x^2 (a+b \text {arcsinh}(c x))^2}{\left (c^2 d x^2+d\right )^2} \, dx\)

\(\Big \downarrow \) 6225

\(\displaystyle \frac {b \int \frac {x (a+b \text {arcsinh}(c x))}{\left (c^2 x^2+1\right )^{3/2}}dx}{c d^2}+\frac {\int \frac {(a+b \text {arcsinh}(c x))^2}{d \left (c^2 x^2+1\right )}dx}{2 c^2 d}-\frac {x (a+b \text {arcsinh}(c x))^2}{2 c^2 d^2 \left (c^2 x^2+1\right )}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {b \int \frac {x (a+b \text {arcsinh}(c x))}{\left (c^2 x^2+1\right )^{3/2}}dx}{c d^2}+\frac {\int \frac {(a+b \text {arcsinh}(c x))^2}{c^2 x^2+1}dx}{2 c^2 d^2}-\frac {x (a+b \text {arcsinh}(c x))^2}{2 c^2 d^2 \left (c^2 x^2+1\right )}\)

\(\Big \downarrow \) 6204

\(\displaystyle \frac {b \int \frac {x (a+b \text {arcsinh}(c x))}{\left (c^2 x^2+1\right )^{3/2}}dx}{c d^2}+\frac {\int \frac {(a+b \text {arcsinh}(c x))^2}{\sqrt {c^2 x^2+1}}d\text {arcsinh}(c x)}{2 c^3 d^2}-\frac {x (a+b \text {arcsinh}(c x))^2}{2 c^2 d^2 \left (c^2 x^2+1\right )}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {\int (a+b \text {arcsinh}(c x))^2 \csc \left (i \text {arcsinh}(c x)+\frac {\pi }{2}\right )d\text {arcsinh}(c x)}{2 c^3 d^2}+\frac {b \int \frac {x (a+b \text {arcsinh}(c x))}{\left (c^2 x^2+1\right )^{3/2}}dx}{c d^2}-\frac {x (a+b \text {arcsinh}(c x))^2}{2 c^2 d^2 \left (c^2 x^2+1\right )}\)

\(\Big \downarrow \) 4668

\(\displaystyle \frac {-2 i b \int (a+b \text {arcsinh}(c x)) \log \left (1-i e^{\text {arcsinh}(c x)}\right )d\text {arcsinh}(c x)+2 i b \int (a+b \text {arcsinh}(c x)) \log \left (1+i e^{\text {arcsinh}(c x)}\right )d\text {arcsinh}(c x)+2 \arctan \left (e^{\text {arcsinh}(c x)}\right ) (a+b \text {arcsinh}(c x))^2}{2 c^3 d^2}+\frac {b \int \frac {x (a+b \text {arcsinh}(c x))}{\left (c^2 x^2+1\right )^{3/2}}dx}{c d^2}-\frac {x (a+b \text {arcsinh}(c x))^2}{2 c^2 d^2 \left (c^2 x^2+1\right )}\)

\(\Big \downarrow \) 3011

\(\displaystyle \frac {2 i b \left (b \int \operatorname {PolyLog}\left (2,-i e^{\text {arcsinh}(c x)}\right )d\text {arcsinh}(c x)-\operatorname {PolyLog}\left (2,-i e^{\text {arcsinh}(c x)}\right ) (a+b \text {arcsinh}(c x))\right )-2 i b \left (b \int \operatorname {PolyLog}\left (2,i e^{\text {arcsinh}(c x)}\right )d\text {arcsinh}(c x)-\operatorname {PolyLog}\left (2,i e^{\text {arcsinh}(c x)}\right ) (a+b \text {arcsinh}(c x))\right )+2 \arctan \left (e^{\text {arcsinh}(c x)}\right ) (a+b \text {arcsinh}(c x))^2}{2 c^3 d^2}+\frac {b \int \frac {x (a+b \text {arcsinh}(c x))}{\left (c^2 x^2+1\right )^{3/2}}dx}{c d^2}-\frac {x (a+b \text {arcsinh}(c x))^2}{2 c^2 d^2 \left (c^2 x^2+1\right )}\)

\(\Big \downarrow \) 2720

\(\displaystyle \frac {2 i b \left (b \int e^{-\text {arcsinh}(c x)} \operatorname {PolyLog}\left (2,-i e^{\text {arcsinh}(c x)}\right )de^{\text {arcsinh}(c x)}-\operatorname {PolyLog}\left (2,-i e^{\text {arcsinh}(c x)}\right ) (a+b \text {arcsinh}(c x))\right )-2 i b \left (b \int e^{-\text {arcsinh}(c x)} \operatorname {PolyLog}\left (2,i e^{\text {arcsinh}(c x)}\right )de^{\text {arcsinh}(c x)}-\operatorname {PolyLog}\left (2,i e^{\text {arcsinh}(c x)}\right ) (a+b \text {arcsinh}(c x))\right )+2 \arctan \left (e^{\text {arcsinh}(c x)}\right ) (a+b \text {arcsinh}(c x))^2}{2 c^3 d^2}+\frac {b \int \frac {x (a+b \text {arcsinh}(c x))}{\left (c^2 x^2+1\right )^{3/2}}dx}{c d^2}-\frac {x (a+b \text {arcsinh}(c x))^2}{2 c^2 d^2 \left (c^2 x^2+1\right )}\)

\(\Big \downarrow \) 6213

\(\displaystyle \frac {2 i b \left (b \int e^{-\text {arcsinh}(c x)} \operatorname {PolyLog}\left (2,-i e^{\text {arcsinh}(c x)}\right )de^{\text {arcsinh}(c x)}-\operatorname {PolyLog}\left (2,-i e^{\text {arcsinh}(c x)}\right ) (a+b \text {arcsinh}(c x))\right )-2 i b \left (b \int e^{-\text {arcsinh}(c x)} \operatorname {PolyLog}\left (2,i e^{\text {arcsinh}(c x)}\right )de^{\text {arcsinh}(c x)}-\operatorname {PolyLog}\left (2,i e^{\text {arcsinh}(c x)}\right ) (a+b \text {arcsinh}(c x))\right )+2 \arctan \left (e^{\text {arcsinh}(c x)}\right ) (a+b \text {arcsinh}(c x))^2}{2 c^3 d^2}+\frac {b \left (\frac {b \int \frac {1}{c^2 x^2+1}dx}{c}-\frac {a+b \text {arcsinh}(c x)}{c^2 \sqrt {c^2 x^2+1}}\right )}{c d^2}-\frac {x (a+b \text {arcsinh}(c x))^2}{2 c^2 d^2 \left (c^2 x^2+1\right )}\)

\(\Big \downarrow \) 216

\(\displaystyle \frac {2 i b \left (b \int e^{-\text {arcsinh}(c x)} \operatorname {PolyLog}\left (2,-i e^{\text {arcsinh}(c x)}\right )de^{\text {arcsinh}(c x)}-\operatorname {PolyLog}\left (2,-i e^{\text {arcsinh}(c x)}\right ) (a+b \text {arcsinh}(c x))\right )-2 i b \left (b \int e^{-\text {arcsinh}(c x)} \operatorname {PolyLog}\left (2,i e^{\text {arcsinh}(c x)}\right )de^{\text {arcsinh}(c x)}-\operatorname {PolyLog}\left (2,i e^{\text {arcsinh}(c x)}\right ) (a+b \text {arcsinh}(c x))\right )+2 \arctan \left (e^{\text {arcsinh}(c x)}\right ) (a+b \text {arcsinh}(c x))^2}{2 c^3 d^2}+\frac {b \left (\frac {b \arctan (c x)}{c^2}-\frac {a+b \text {arcsinh}(c x)}{c^2 \sqrt {c^2 x^2+1}}\right )}{c d^2}-\frac {x (a+b \text {arcsinh}(c x))^2}{2 c^2 d^2 \left (c^2 x^2+1\right )}\)

\(\Big \downarrow \) 7143

\(\displaystyle \frac {2 \arctan \left (e^{\text {arcsinh}(c x)}\right ) (a+b \text {arcsinh}(c x))^2+2 i b \left (b \operatorname {PolyLog}\left (3,-i e^{\text {arcsinh}(c x)}\right )-\operatorname {PolyLog}\left (2,-i e^{\text {arcsinh}(c x)}\right ) (a+b \text {arcsinh}(c x))\right )-2 i b \left (b \operatorname {PolyLog}\left (3,i e^{\text {arcsinh}(c x)}\right )-\operatorname {PolyLog}\left (2,i e^{\text {arcsinh}(c x)}\right ) (a+b \text {arcsinh}(c x))\right )}{2 c^3 d^2}+\frac {b \left (\frac {b \arctan (c x)}{c^2}-\frac {a+b \text {arcsinh}(c x)}{c^2 \sqrt {c^2 x^2+1}}\right )}{c d^2}-\frac {x (a+b \text {arcsinh}(c x))^2}{2 c^2 d^2 \left (c^2 x^2+1\right )}\)

input
Int[(x^2*(a + b*ArcSinh[c*x])^2)/(d + c^2*d*x^2)^2,x]
 
output
-1/2*(x*(a + b*ArcSinh[c*x])^2)/(c^2*d^2*(1 + c^2*x^2)) + (b*(-((a + b*Arc 
Sinh[c*x])/(c^2*Sqrt[1 + c^2*x^2])) + (b*ArcTan[c*x])/c^2))/(c*d^2) + (2*( 
a + b*ArcSinh[c*x])^2*ArcTan[E^ArcSinh[c*x]] + (2*I)*b*(-((a + b*ArcSinh[c 
*x])*PolyLog[2, (-I)*E^ArcSinh[c*x]]) + b*PolyLog[3, (-I)*E^ArcSinh[c*x]]) 
 - (2*I)*b*(-((a + b*ArcSinh[c*x])*PolyLog[2, I*E^ArcSinh[c*x]]) + b*PolyL 
og[3, I*E^ArcSinh[c*x]]))/(2*c^3*d^2)
 

3.3.36.3.1 Defintions of rubi rules used

rule 27
Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

rule 216
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[b, 2]))*A 
rcTan[Rt[b, 2]*(x/Rt[a, 2])], x] /; FreeQ[{a, b}, x] && PosQ[a/b] && (GtQ[a 
, 0] || GtQ[b, 0])
 

rule 2720
Int[u_, x_Symbol] :> With[{v = FunctionOfExponential[u, x]}, Simp[v/D[v, x] 
   Subst[Int[FunctionOfExponentialFunction[u, x]/x, x], x, v], x]] /; Funct 
ionOfExponentialQ[u, x] &&  !MatchQ[u, (w_)*((a_.)*(v_)^(n_))^(m_) /; FreeQ 
[{a, m, n}, x] && IntegerQ[m*n]] &&  !MatchQ[u, E^((c_.)*((a_.) + (b_.)*x)) 
*(F_)[v_] /; FreeQ[{a, b, c}, x] && InverseFunctionQ[F[x]]]
 

rule 3011
Int[Log[1 + (e_.)*((F_)^((c_.)*((a_.) + (b_.)*(x_))))^(n_.)]*((f_.) + (g_.) 
*(x_))^(m_.), x_Symbol] :> Simp[(-(f + g*x)^m)*(PolyLog[2, (-e)*(F^(c*(a + 
b*x)))^n]/(b*c*n*Log[F])), x] + Simp[g*(m/(b*c*n*Log[F]))   Int[(f + g*x)^( 
m - 1)*PolyLog[2, (-e)*(F^(c*(a + b*x)))^n], x], x] /; FreeQ[{F, a, b, c, e 
, f, g, n}, x] && GtQ[m, 0]
 

rule 3042
Int[u_, x_Symbol] :> Int[DeactivateTrig[u, x], x] /; FunctionOfTrigOfLinear 
Q[u, x]
 

rule 4668
Int[csc[(e_.) + Pi*(k_.) + (Complex[0, fz_])*(f_.)*(x_)]*((c_.) + (d_.)*(x_ 
))^(m_.), x_Symbol] :> Simp[-2*(c + d*x)^m*(ArcTanh[E^((-I)*e + f*fz*x)/E^( 
I*k*Pi)]/(f*fz*I)), x] + (-Simp[d*(m/(f*fz*I))   Int[(c + d*x)^(m - 1)*Log[ 
1 - E^((-I)*e + f*fz*x)/E^(I*k*Pi)], x], x] + Simp[d*(m/(f*fz*I))   Int[(c 
+ d*x)^(m - 1)*Log[1 + E^((-I)*e + f*fz*x)/E^(I*k*Pi)], x], x]) /; FreeQ[{c 
, d, e, f, fz}, x] && IntegerQ[2*k] && IGtQ[m, 0]
 

rule 6204
Int[((a_.) + ArcSinh[(c_.)*(x_)]*(b_.))^(n_.)/((d_) + (e_.)*(x_)^2), x_Symb 
ol] :> Simp[1/(c*d)   Subst[Int[(a + b*x)^n*Sech[x], x], x, ArcSinh[c*x]], 
x] /; FreeQ[{a, b, c, d, e}, x] && EqQ[e, c^2*d] && IGtQ[n, 0]
 

rule 6213
Int[((a_.) + ArcSinh[(c_.)*(x_)]*(b_.))^(n_.)*(x_)*((d_) + (e_.)*(x_)^2)^(p 
_.), x_Symbol] :> Simp[(d + e*x^2)^(p + 1)*((a + b*ArcSinh[c*x])^n/(2*e*(p 
+ 1))), x] - Simp[b*(n/(2*c*(p + 1)))*Simp[(d + e*x^2)^p/(1 + c^2*x^2)^p] 
 Int[(1 + c^2*x^2)^(p + 1/2)*(a + b*ArcSinh[c*x])^(n - 1), x], x] /; FreeQ[ 
{a, b, c, d, e, p}, x] && EqQ[e, c^2*d] && GtQ[n, 0] && NeQ[p, -1]
 

rule 6225
Int[((a_.) + ArcSinh[(c_.)*(x_)]*(b_.))^(n_.)*((f_.)*(x_))^(m_)*((d_) + (e_ 
.)*(x_)^2)^(p_), x_Symbol] :> Simp[f*(f*x)^(m - 1)*(d + e*x^2)^(p + 1)*((a 
+ b*ArcSinh[c*x])^n/(2*e*(p + 1))), x] + (-Simp[f^2*((m - 1)/(2*e*(p + 1))) 
   Int[(f*x)^(m - 2)*(d + e*x^2)^(p + 1)*(a + b*ArcSinh[c*x])^n, x], x] - S 
imp[b*f*(n/(2*c*(p + 1)))*Simp[(d + e*x^2)^p/(1 + c^2*x^2)^p]   Int[(f*x)^( 
m - 1)*(1 + c^2*x^2)^(p + 1/2)*(a + b*ArcSinh[c*x])^(n - 1), x], x]) /; Fre 
eQ[{a, b, c, d, e, f}, x] && EqQ[e, c^2*d] && GtQ[n, 0] && LtQ[p, -1] && IG 
tQ[m, 1]
 

rule 7143
Int[PolyLog[n_, (c_.)*((a_.) + (b_.)*(x_))^(p_.)]/((d_.) + (e_.)*(x_)), x_S 
ymbol] :> Simp[PolyLog[n + 1, c*(a + b*x)^p]/(e*p), x] /; FreeQ[{a, b, c, d 
, e, n, p}, x] && EqQ[b*d, a*e]
 
3.3.36.4 Maple [F]

\[\int \frac {x^{2} \left (a +b \,\operatorname {arcsinh}\left (c x \right )\right )^{2}}{\left (c^{2} d \,x^{2}+d \right )^{2}}d x\]

input
int(x^2*(a+b*arcsinh(c*x))^2/(c^2*d*x^2+d)^2,x)
 
output
int(x^2*(a+b*arcsinh(c*x))^2/(c^2*d*x^2+d)^2,x)
 
3.3.36.5 Fricas [F]

\[ \int \frac {x^2 (a+b \text {arcsinh}(c x))^2}{\left (d+c^2 d x^2\right )^2} \, dx=\int { \frac {{\left (b \operatorname {arsinh}\left (c x\right ) + a\right )}^{2} x^{2}}{{\left (c^{2} d x^{2} + d\right )}^{2}} \,d x } \]

input
integrate(x^2*(a+b*arcsinh(c*x))^2/(c^2*d*x^2+d)^2,x, algorithm="fricas")
 
output
integral((b^2*x^2*arcsinh(c*x)^2 + 2*a*b*x^2*arcsinh(c*x) + a^2*x^2)/(c^4* 
d^2*x^4 + 2*c^2*d^2*x^2 + d^2), x)
 
3.3.36.6 Sympy [F]

\[ \int \frac {x^2 (a+b \text {arcsinh}(c x))^2}{\left (d+c^2 d x^2\right )^2} \, dx=\frac {\int \frac {a^{2} x^{2}}{c^{4} x^{4} + 2 c^{2} x^{2} + 1}\, dx + \int \frac {b^{2} x^{2} \operatorname {asinh}^{2}{\left (c x \right )}}{c^{4} x^{4} + 2 c^{2} x^{2} + 1}\, dx + \int \frac {2 a b x^{2} \operatorname {asinh}{\left (c x \right )}}{c^{4} x^{4} + 2 c^{2} x^{2} + 1}\, dx}{d^{2}} \]

input
integrate(x**2*(a+b*asinh(c*x))**2/(c**2*d*x**2+d)**2,x)
 
output
(Integral(a**2*x**2/(c**4*x**4 + 2*c**2*x**2 + 1), x) + Integral(b**2*x**2 
*asinh(c*x)**2/(c**4*x**4 + 2*c**2*x**2 + 1), x) + Integral(2*a*b*x**2*asi 
nh(c*x)/(c**4*x**4 + 2*c**2*x**2 + 1), x))/d**2
 
3.3.36.7 Maxima [F]

\[ \int \frac {x^2 (a+b \text {arcsinh}(c x))^2}{\left (d+c^2 d x^2\right )^2} \, dx=\int { \frac {{\left (b \operatorname {arsinh}\left (c x\right ) + a\right )}^{2} x^{2}}{{\left (c^{2} d x^{2} + d\right )}^{2}} \,d x } \]

input
integrate(x^2*(a+b*arcsinh(c*x))^2/(c^2*d*x^2+d)^2,x, algorithm="maxima")
 
output
-1/2*a^2*(x/(c^4*d^2*x^2 + c^2*d^2) - arctan(c*x)/(c^3*d^2)) + integrate(b 
^2*x^2*log(c*x + sqrt(c^2*x^2 + 1))^2/(c^4*d^2*x^4 + 2*c^2*d^2*x^2 + d^2) 
+ 2*a*b*x^2*log(c*x + sqrt(c^2*x^2 + 1))/(c^4*d^2*x^4 + 2*c^2*d^2*x^2 + d^ 
2), x)
 
3.3.36.8 Giac [F]

\[ \int \frac {x^2 (a+b \text {arcsinh}(c x))^2}{\left (d+c^2 d x^2\right )^2} \, dx=\int { \frac {{\left (b \operatorname {arsinh}\left (c x\right ) + a\right )}^{2} x^{2}}{{\left (c^{2} d x^{2} + d\right )}^{2}} \,d x } \]

input
integrate(x^2*(a+b*arcsinh(c*x))^2/(c^2*d*x^2+d)^2,x, algorithm="giac")
 
output
integrate((b*arcsinh(c*x) + a)^2*x^2/(c^2*d*x^2 + d)^2, x)
 
3.3.36.9 Mupad [F(-1)]

Timed out. \[ \int \frac {x^2 (a+b \text {arcsinh}(c x))^2}{\left (d+c^2 d x^2\right )^2} \, dx=\int \frac {x^2\,{\left (a+b\,\mathrm {asinh}\left (c\,x\right )\right )}^2}{{\left (d\,c^2\,x^2+d\right )}^2} \,d x \]

input
int((x^2*(a + b*asinh(c*x))^2)/(d + c^2*d*x^2)^2,x)
 
output
int((x^2*(a + b*asinh(c*x))^2)/(d + c^2*d*x^2)^2, x)